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Abstract 



If^ We further clarify how scalar metric perturbations are amplified in an 

On infiationary cosmology. We first construct a simple, analytic model of an in- 

fiationary cosmology in which the expansion scale factor evolves continuously 



Oh 



from an infiationary era to a radiation-dominated era. From this model, it 



is clear to see how scalar perturbations are amplified. Second, we examine 
the recent claims of Grishchuk, and the reply by Deruelle and Mukhanov, 
regarding the evolution of scalar perturbations through an abrupt transition 
in the equation of state of the cosmological fiuid. We demonstrate that the 
"standard results" regarding the amplification of scalar, density perturbations 
from infiation are valid. 

o 

On 
O 

XT) I. INTRODUCTION 

^ , The goal of this paper is to understand the recent claims by Grishchuk [1-3] that the stan- 

^ dard results on the amplification of scalar metric perturbations in an infiationary cosmology 

(3X) are incorrect. By standard result, we mean that the magnitude of the scalar perturbations 
is amplified by a factor determined by the change in the equation of state w = p/ p of the 
cosmological fiuid, where p and p are the background pressure and energy density. Crudely, 
the amplification is proportional to 

1 + ^ (1.1) 

when the perturbations evolve through consecutive expansion eras in which the equation 
of state of the cosmological fiuid is given by w^. For |1 + W-\ ^ 1, as occurs in an 
infiationary epoch, this amplification may be large. 

It was argued in references [1-3] that the amplification of density perturbations was much 
weaker than commonly assumed. However, this claim is incorrect, as has been pointed out 
by Deruelle and Mukhanov [4]. There it was shown that the error in Grishchuk's argument 
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was due to a subtlety in the junction conditions necessary to evolve the density perturbations 
across an abrupt transition in the equation of state of the cosmological fluid at the end of 
the inflationary epoch. 

While the issue would seem to be settled, we found it useful to study a cosmologi- 
cal scenario which made no recourse to junction conditions. Hence, we constructed a toy 
cosmological model which evolves smoothly from an inflationary- to radiation-dominated ex- 
pansion era. Although the evolution of density perturbations in this model must be carried 
out numerically, implementation of the code is a straight forward task. Hence, we confirmed 
that the standard results (as presented in [5] for example) were indeed valid. We hope 
that this paper will contribute towards the understanding of the mechanism by which scalar 
metric perturbations are amplified in an infiationary cosmology. 

The contents of this paper are as follows. In section H, we construct a toy cosmological 
model and numerically evolve the scalar metric perturbations across the infiation-radiation 
transition. In section III, we point out the shortcomings in Grishchuk's argument, making 
reference to the work of Deruelle and Mukhanov. In section IV we turn to the cosmological 
significance of scalar metric perturbations. We demonstrate that scalar perturbations do 
not vanish in the de Sitter limit, contrary to Grishchuk's claim, but in fact are strongly 
amplified. We conclude in section V. Throughout this paper we use units in which the 
speed of light c = I. 



II. TOY MODEL INFLATIONARY COSMOLOGY 

In this section we construct an analytic, toy model infiationary cosmology. This model 
possesses the generic features of an infiationary cosmology; an arbitrarily long vacuum- 
dominated expansion era followed by a smooth transition to a radiation-dominated era, and 
strong amplification of scalar perturbations. We refer to this model as a "toy" because we 
have not made contact with a sensible particle physics model; this model is proposed for 
simple, demonstrative purposes. 



A. Cosmological Model 

Consider a spatially-fiat, isotropic FRW space-time with line element 

ds^ = -dt^ + a^{t)[dx^ + dy^ + dz\ (2.1) 

Here, t is the physical time, and ait) is the expansion scale factor. The matter content of 
the space-time is given by 

3^-2 I 

'^'^ = 8^^K^- ^^-^^ 

Here, ai is a dimensionless constant, and L is a length scale. Examining the above equation, 
we see that for small values of the scale factor, a ^ ai, the energy density is approximately 
constant, as in a de Sitter, vacuum-dominated space-time. For large values of the scale- 
factor, a ^ ai, the energy density decays as a""*, as in a radiation-dominated space-time. 
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Thus, we have a "two-fluid" cosmology which evolves from vacuum- to radiation-domination 
with increasing a. We may adjust the energy density around the time of the transition from 
vacuum- to radiation-domination, by varying ai. A simple extension of this toy model is to 
add another component which decays as a~"^, as in a dust-dominated space-time. In such a 
case, the energy density would be 



where 02 is a dimensionless constant which determines the time of the transition from 
radiation- to dust-domination. For simplicity, we will only consider the toy model energy 
density given by equation (2.2). 

The remaining details of the cosmology may be obtained from the Einstein equations. 
The Hubble rate is given by 

H = ^^ = L-'{a\ + a*)-'l\ (2.3) 
a at 

The pressure of the cosmological fluid is 

a'^ — 'ia\ la'* — 'ia\ 
STTG{a\ + a^y " 3 a4 + «^ ^' 

From the equation for the energy density and pressure, we see that the inequality /9-|-3p < 
is satisfled for a < ai. Thus, the space-time will inflate for a < ai. We are also interested in 
the quantities 

p l/a4-3at 
w = — = -[ — j- 

p 3^ a'^ -\- aj 

_ dp dpi dt 1 

7 = -Uro,y = ^ = -- + 2W. (2.5) 

Here, w gives the equation of state of the cosmological fluid, and 7 is the adiabatic index. 

An advantage of this toy model, for instructional purposes, is that it makes no refer- 
ence to a scalar held or potential which may compose the vacuum-dominated cosmological 
fluid. As well, all quantities are expressed in terms of the expansion scale factor a and the 
dimensionless constant ai. There is no need to obtain the time dependence of the scale 
factor. 



B. Adiabatic Density Perturbations 

We consider scalar perturbations to the spatially flat FRW metric [5] 

ds"^ = a^{ri)[-dri^{l + 2</)) + 2B,,dx'dr] + dx'dx'{2E,,, + (1 - 2^')%)], (2.6) 

where t] is the conformal time. In the above line element, the derivatives are with respect 
to the spatially flat background 3— metric. We may construct combinations of the scalar 
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quantities (/>, B which are invariant under a general coordinate transformation x" — > 

x" = x" -\- ^" . These invariant combinations 

$ = (j) + -WB - E')a']' 
a 

m = 4^--{B - E') (2.7) 
a 

help simplify the evolution equations for density perturbations. Here ' = d/drj. Similarly, we 
may consider the perturbed Einstein equations under a general coordinate transformation 
and construct gauge invariant quantities. These equations may be written in terms of $ and 
^, and gauge invariant combinations of the source stress-energy tensor. For a symmetric 
source stress-energy tensor, a further simplification occurs, as the field equations require 
$ = ^. The quantity $ = is referred to as the gauge invariant potential, and 

characterizes the amplitude of density perturbations. On scales below the Hubble radius, $ 
plays a role similar to the Newtonian potential. 

The evolution of the gauge invariant potential for the case of adiabatic perturbations is 
given by the equation 

^ + 37/(1 + 7)— - 7V2$ + [2— + (1 + 37)7/2]$ = o. (2.8) 

By adiabatic perturbations, we mean perturbations to the cosmological fiuid 5p and 5p such 
that 5p = ^Sp. In the above evolution equation, rj is the conformal expansion time, and 
H = dlna/drj. A clear derivation of this evolution equation is found in the review by 
Mukhanov, Feldman, and Brandenberger [5]. We would like to follow the evolution of $ 
from the vacuum- to radiation-dominated era. 

We will recast the evolution equation using the expressions dt = adrj and da/dt = aH to 
relate the conformal and physical time to the scale factor. We obtain 

+ -(tt + 37 - t:^)-, JTT^V^^ + —3(7 - w)<!> = 0. (2.9) 



da"" a'2 ' 2 ' da a 

Thus, we have eliminated the need for a time coordinate, expressing the evolution in terms 
of the expansion scale factor. Note that this form of the evolution equation is generally 
valid; we have not used any results specific to our toy cosmology in arriving at (2.9). 
We next decompose the gravitational potential into spatial Fourier harmonics. 

$(a,f) = j d^M{a,k)e'''-^ (2.10) 

Thus we replace $ — > $ and — > —k^ in equation (2.9). This decomposition will allow 
us to follow the evolution of a single mode through the vacuum- and radiation-dominated 
expansion eras. 

Let us qualitatively examine the evolution equation (2.9). At early and late times, we 
may use equations (2.5) to simplify the form of the evolution equation. 



For a«ai, - - {-^k' + -)<^ = 

aa^ a da 3a^ 
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At early times, 7 < t/;, and $ has a rapidly growing mode. Physically, the quickly changing 
negative vacuum pressure of the space-time couples to the gauge invariant potential, serving 
as a strong anti-damping or "pump" source. Hence, we expect rapid growth in $ at early 
times. This is the mechanism by which the scalar perturbations are amplified. The equation 
of state is ~ |(a/ai)'* ^ 1, indicating strong amplification according to the expression 

(1.1). At late times, 7 ~ t/;, and the coefficient of the term linear in $ is negligibly small for 
k^L^ <^ a\l , indicating a constant solution. Hence, we expect the rapid growth to cease 
and the amplitude of $ to stabilize at late times. While this is indeed what happens, we are 
not interested in a crude, qualitative estimate. We proceed to the next section for an exact 
solution to the evolution. 



C. Numerical Solution 



The evolution of the gauge invariant potential as a function of the expansion scale factor 
may easily be obtained numerically. It is a straight forward exercise to convert equation (2.9) 
into a pair of first order differential equations. Specifying the initial data ($, c/ln $/c/ln a) 
at some initial time a^, we may integrate to obtain $ at a later time aj. The results are 
presented in figure 1. Curves showing the evolution of $ are presented in figure 1, for 
initial values of c/ln$/c/lna = (0, 1, 10), with k <^ representing scalar fiuctuations with 
wavelengths longer than the Hubble length. The potential grows rapidly, as $ oc a"* in the 
vacuum-dominated era. However, upon the onset of the radiation-dominated era, for a ^ ai, 
the growth stops and $ remains constant. The evolution is just as we expected from our 
qualitative discussion above. Hence, we have constructed a toy cosmological model in which 
the magnitude of scalar perturbations is strongly amplified. 



III. STANDARD TECHNIQUES 



In this section we examine the recent claims of Grishchuk [1-3], and the reply by Deruelle 
and Mukhanov [4], regarding the evolution of scalar perturbations across a phase transition. 
We demonstrate that the "standard results" regarding the amplification of scalar, density 
perturbations from infiation are valid. 



A. as a Conserved Quantity 



A standard technique for estimating the amplitude of the gauge invariant potential $ is 
to use the integral of the equation of motion, (. As has been well demonstrated elsewhere 
[5], the quantity 

6 I + W 
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is approximately constant in time for modes k <^ H. Thus, one may use C as a sort of 
conserved quantity. Applying this technique to the case of our toy model cosmology, we 
obtain 

= ^ I = {.(°nM±4+ ;^+ "■;''■'*■ ). (3.2) 

Here, we have assumed that c/ln $/c/ln a ^ 1. A comparison of the evolution of $ using 
equation (3.2) with the exact, numerical result is presented in figure 2. We see that the use 
of C as a conserved quantity to obtain the evolution of $ is well justified. 

It has been claimed by Grishchuk [2] that C = foi" times, and therefore conveys no 
information as a conserved quantity. We will examine this argument closely. In equation 
(38') of [2], the evolution equation for the synchronous gauge scalar quantity /i is given as 



)Un-^ = X. (3.3) 



Here 7g = 1 — H' /H^, n is the mode number (referred to in these notes as k), and 

/i = -i-^(V'' + 'H7GV') (3.4) 



with defined by equation (2.6). Next, taking the definition of ( given by [5], expressed in 
terms of synchronous gauge variables, Grishchuk writes 

CMFB = 7r^— — \a -fG[ — ^) • 3.5 



The argument proceeds as follows. Take the long wavelength limit, for which C, is claimed 
to be conserved, letting n — > in equation (3.3). Insert the resulting expression for X into 
the equation defining C, to obtain 

Cmfb = TT-;- (3.6) 

2n^ 

Note that (3.3) is simply the equation of motion for the scalar perturbations, with X = 0. 
Hence, Grishchuk claims, (mfb = 0. However, the error in this argument is apparent. The 
n — > limit has not been taken consistently. Returning to (3.3), using X = 0, we substitute 
the expression for n^/i/a^/7G into (3.5) to obtain 

Cmfb = ^. — • (3.7) 



Relating [i to our more familiar gauge invariant potential $, we find that this equation is 
equivalent to (3.1). Hence, this exercise has simply used the equations of motion to recast 
the form of C,. 

The quantity C, may be zero for specially chosen data. For a given background space- 
time, specifying the matter content or the expansion scale factor, one may always choose 
a set of initial data ($i, c/ln $i/c/ln a^) satisfying the equation of motion (2.9), for which 
C = 0. However, the equations of motion do not demand C = 0. 
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B. Matching Conditions for the Evolution Across an Abrupt Transition 



In the recent work of Deruelle and Mukhanov [4], the junction conditions were obtained 
which are required to patch together the evolution of the gauge invariant potential across 
a phase transition in which the equation of state of the cosmological fluid changes discon- 
tinuously. The scenario considered is the following. The gauge invariant potential evolves 
according to equation (2.9), until the adiabatic index 7 and the equation of state w change 
discontinuously, after which the evolution is again described by (2.9). To match the solutions 
$_ and $_|_ (the it refers the value on either side of the transition) across the transition, 
they required the continuity of the total energy and momentum of the cosmological fluid. 
Defining the transition to occur when the total pressure, the background plus perturbations, 
reaches a critical value, they showed that the surface of transition is not in general a surface 
of constant time. Transforming into a coordinate system in which the transition surface oc- 
curs at constant time, they obtain the matching conditions, in terms of synchronous gauge 
quantities ?/> and E: 

In the synchronous gauge, chosen here to make contact with Grishchuk's work, (j) = B = 
in equations (2.6) and (2.7)). We see that the matching conditions reproduce the result that 
( serves as a conserved quantity for long wavelength perturbations, and is constant across 
the transition. 

As an application of these matching conditions, consider a power-law expansion epoch, 
such that ^ = w. Given initial conditions ($i, c/ln $i/c/ln a^), the evolution for $ is 

for long wavelength modes. Employing the matching conditions (3.8) to continue the evo- 
lution into another power-law expansion epoch we find 

c/ln$_i_ /Wa_ — W-\ -|- 1 \ c/ln $_ 

+ {^^)^^ • (3-10) 



din a ^ -\- 1 ' ^W-^\' din a 

For |1 -|- W- \ ^ 1, as occurs at the end of an infiationary epoch, the derivative d^ j da term 
will receive a substantial amplification. The evolution of $ through vacuum- {w = —\ ^ 
< ^ 1) and radiation-dominated {w = 1/3) expansion epochs is displayed in figure 3. 
While no substantial growth occurs for the long wavelength modes during either epoch, the 
abrupt transition leads to a strong amplification. The late-time amplification obtained by 
modeling the transition from vacuum- to radiation-dominated expansion by a discontinuous 
jump in the equation of state of the cosmological fiuid agrees with the exact result obtained 
for the cosmological model used in section II. 

Grishchuk has proposed a different set of matching conditions for the evolution of the 
gauge invariant potential across an abrupt transition. In terms of Grishchuk's quantities, 
Deruelle and Mukhanov's matching conditions (3.8) are written as 
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[-h + n-2m;g]± = o [c - ^^7^77— -v^i-/. + n-^m\Q)]± = 0, (3.11) 

Where Q{x) are the scalar harmonics on W. In terms of our notation, h = —tp/2 and 
hi = n'^E/2. However, Grishchuk [1] has matched the scalar quantities /i, hi on a surface 
of constant time which happens not to be a surface of constant energy. This amounts 
to requiring that the scalar quantities /i, hi and their time derivatives match across the 
transition: 

[h]± = [h']± = N± = [h'i]± = 0. (3.12) 

These conditions are clearly inequivalent to those obtained by Deruelle and Mukhanov, 
which return the standard result. 

An additional set of matching conditions are discussed by Grishchuk. Following equation 
(48) of [1], the evolution of the scalar perturbations across the abrupt transition requires 

Mi = H± = (3.13) 
as the matching conditions. The quantity /i has been defined in equation (3.4), and 

v = -n^h + h'i = + E'). (3.14) 

We see that matching v is consistent with (3.12). Next, using (3.7) for ( in terms of /i, and 
(3.4) for /i in terms of /i, /i', we rewrite the additional matching condition as 

M± = [-2av^C]± = [t7^(/^' + n^Gh)U = 0. (3.15) 



It is clear from the middle term above that this condition is inequivalent to the correct 
condition (3.8), since ( is conserved for long wavelengths, whereas 7g changes discontinuously 
across the abrupt transition. From the term on the right above, we see that matching 
/i is inequivalent to matching /i, h' . Hence, the conditions (3.13) are inequivalent to the 
conditions actually employed by Grishchuk, (3.12), and the correct conditions (3.8). 

The evolution of $ through consecutive power-law expansion epochs, using both condi- 
tions (3.12) and (3.13), are shown in figure 3. Clearly, these choices of matching conditions 
are incorrect, as they disagree with both the exact results, and those results obtained through 
other techniques. 

We take this opportunity to comment briefiy on the matching conditions for cosmological 
perturbations at a transition in the equation of state prescribed by Ratra [7]. The condition 
given for the evolution of the density perturbation 5p^ equations (3.10) and (4.8) in [7] agree 
with the results of Deruelle and Mukhanov [4] only in the longitudinal gauge. The results do 
not agree, and are incorrect when applied in the synchronous gauge. The reason is that the 
density perturbation itself is a gauge dependent quantity; Ratra's condition for the joining 
of 5p IS not gauge invariant. 
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IV. STANDARD RESULTS 



In this brief section we make contact with the standard results on the amphtude of per- 
turbations produced by inflation. We first discuss how the amplitude of scalar perturbations 
is determined. Second, we consider the amplitude of perturbations in the limiting case of a 
de Sitter infiationary era. 



A. Observable Quantities 

The physical significance of the gauge invariant potential may be understood if we con- 
sider the temperature anisotropy induced in the cosmic microwave background (CMB). On 
large scales, the dominant source of the temperature anisotropy due to the gravitational 
potential at the surface of last scattering [6] is 

ST 1 

— (0) ^ -$(?7i„|?7o-??is|0). (4.1) 

Thus, the CMB anisotropy pattern is determined by the gauge invariant potential $. In 
this equation the potential is evaluated at the time of last scattering rjis^ where the present 
time is r/o, Vo ~ Vis is the conformal distance to the surface of last scattering, and is a unit 
vector on the sphere. 

From our experience in the previous sections, we may use the quantity ( to evolve the 
potential $ from the infiationary to matter-dominated era. Hence, for long wavelengths, 
in the matter era is given by 

$™ ^ r ^ (4-2) 

5 1 + w, 

where Wi^ $i are evaluated in the infiationary era. Next we must set the initial amplitude 
of the potential. 

In a realistic model of infiation, scalar metric perturbations are generated according to 
the following crude scheme. The vacuum energy of a homogeneous scalar field (j)o drives the 
infiationary expansion; quantum mechanical fiuctuations 5(f around the background scalar 
field produce scalar metric fiuctuations $ [5,6]; the scalar perturbations are amplified by the 
vacuum-dominated expansion. As we have already examined the subsequent evolution of 
the scalar perturbations, we need to fix the initial conditions for the potential $. We refer 
to Einstein's equations to express the fiuctuations around the background source in terms 
of metric perturbations (see [5]), 

m + = 47TG<P',5ip ^ <^m^l'H^\. (4.3) 

5 00 

where we have used the relationship 1 -\- w = (SttG /3)'H~'^ (j)^ for the scalar field-dominated 
equation of state. The amplitude of 5(f is determined by the equal time, two-point function 
{5(f[x)5(f[x')) in the appropriate quantum mechanical vacuum state. Hence, the amplitude 
of the CMB anisotropy, and the gauge invariant potential, is determined by the amplitude 
of the quantum mechanical fiuctuations in the background scalar field. Equation (4.3) for 
$m agrees with the standard result for the amplitude of perturbations from infiation [5,8]. 
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B. The Limit of de Sitter Inflation 



We now examine the case of a de Sitter inflationary era, taken as the hmit of a general 
power-law expansion era. Consider a power-law expansion scale factor a{ri) oc ri~^~^ where 
r] is the conformal time, as — > 0. In this well-studied scenario [6], the time variation of the 
background scalar held 

vanishes in the de Sitter, 5^0 limit. Returning to (4.3), the term which flxes the amplitude 
of the potential during the de Sitter era, the product </>o^'/?, appears to vanish. However, let 
us also examine the rms amplitude of 5(f. In this scenario, 

(V) = r dkk'\H^ll{kr^)\\ (4.5) 

In the de Sitter limit, it is well known that the rms amplitude diverges logarithmically, due 
to the behavior of the integrand at small k. Examining the integral more closely, we flnd 
that in the limit — > 0, the product (/)q(^(/7^)^/^ does not vanish. Hence, the term in (4.3) 
which flxes the amplitude of the potential does not vanish in the de Sitter limit. 

We may now return to the claims by Grishchuk [1-3] that the magnitude of scalar 
perturbations vanishes in the limit of a de Sitter inflationary era. The fault in the argument 
is that the divergence in the rms amplitude of the fluctuations in the scalar held has been 
ignored. As we have demonstrated above, the amplitude of scalar perturbations is in fact 
flnite during the de Sitter epoch. However, these scalar perturbations experience inflnite 
ampliflcation at the close of the inflationary era, since the factor 1 -|- in (1.1) vanishes. 
Note that the calculation of the evolution of scalar perturbations has been carried out 
within the context of linearized gravity, under the assumption that all perturbations are 
small, |$| ^ 1. Hence, there is a limit to how closely the inflationary era may consistently 
approximate a de Sitter space-time. Generically, this restriction is satisfled by inflationary 
cosmologies for which the amplitude of the induced temperature anisotropics is consistent 
with observations. 

V. SUMMARY 

We have conflrmed the validity of the standard results on the evolution of scalar pertur- 
bations in an inflationary cosmology. Our main result is the demonstration that the gauge 
invariant potential $ is strongly amplifled during the course of evolution from a vacuum- to 
radiation-dominated expansion era. Next, we have demonstrated that C, the integral of the 
time dependent portion of the equation of motion, may be used to excellent approximation 
as a tool for computing the evolution of long wavelength modes of $. 

We have also discussed the recent claims of Grishchuk [1-3] that scalar perturbations 
are not strongly amplifled in an inflationary cosmology. First, we demonstrated that the 
claim C = for all times is incorrect. Second, making reference to the work of Deruelle and 
Mukhanov [4], we have demonstrated the flaw in Grishchuk's matching conditions for the 
evolution of the gauge invariant potential across an abrupt transition. The standard result 
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is consistent with the conservation of ( for long wavelength perturbations. Third, we have 
shown that scalar perturbations do not vanish in the limit of a de Sitter inflationary epoch. 
Rather, the perturbations are strongly amplified as the limit is approached, as given in the 
standard literature [4-6]. 
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FIGURES 
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FIG. 1. The evolution of the gauge invariant potential ^ in the toy cosmological model con- 
structed in these notes is displayed. The evolution is computed exactly, for three sets of initial 
conditions, for long wavelength modes. The expansion is vacuum-dominated for a < ai, and 
radiation-dominated for ai < a. 
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FIG. 2. The evolution of the gauge invariant potential <i> in the toy cosmological model con- 
structed in these notes is displayed. The solid curve shows the exact evolution. The dashed curve 
shows the evolution obtained using as a conserved quantity. The disagreement for a < ai occurs 
because we have ignored the contribution of the derivative term rfln <i>/(iln a in (. The agreement 
is excellent at late times, when |(iln <i>/(iln a| <^ 1 and the derivative term is negligible. 
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FIG. 3. The evolution of the gauge invariant potential <i> in a cosmological model constructed 
from consecutive, power-law expansion epochs, matching the potential across the abrupt transition. 
The evolution is computed exactly, for three sets of initial conditions, for long wavelength modes. 
The matching conditions obtained by Deruelle and Mukhanov, using as a conserved quantity, 
were used to evolve <i> across the transition. The light, solid lines show the evolution of <i> using 
the matching conditions (3.12) used by Grishchuk. The light, dotted lines show the evolution 
of <i> using the matching conditions (3.13) discussed by Grishchuk. The light, long-short dashed 
lines show the exact evolution of <i> from figure 1. The expansion is vacuum-dominated for a < ai 
(w = — l-|-|(aj/ai)'*), and radiation-dominated for ai < a (w = 1/3). Similar results were obtained 
for a toy model which evolves from vacuum- to dust-dominated expansion. 
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